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ABSTRACT 

In order to reproduce the low mass end of the steUar mass function, most cur- 
rent models of galaxy evolution invoke very efficient supernova feedback. This solution 
seems to suffer from several shortcomings however, like predicting too little star for- 
mation in low mass galaxies at z=0. In this work, we explore modifications to the star 
formation (SF) law as an alternative solution to achieve a match to the stellar mass 
function. This is done by applying semi- analytic models based on De Lucia & Blaizot, 
but with varying SF laws, to the Millennium and Millennium-II simulations, within 
the formalism developed by Neistein & Weinmann. Our best model includes lower 
SF efficiencies than predicted by the Kennicutt- Schmidt law at low stellar masses, no 
sharp threshold of cold gas mass for SF, and a SF law that is independent of cosmic 
time. These simple modifications result in a model that is more successful than cur- 
rent standard models in reproducing various properties of galaxies less massive than 
lO^^Mo. The improvements include a good match to the observed auto-correlation 
function of galaxies, an evolution of the stellar mass function from z = 3 to z = 
similar to observations, and a better agreement with observed specific star formation 
rates. However, our modifications also lead to a dramatic overprediction of the cold 
mass content of galaxies. This shows that finding a successful model may require fine- 
tuning of both star formation and supernovae feedback, as well as improvements on 
gas cooling, or perhaps the inclusion of a yet unknown process which efficiently heats 
or expels gas at high redshifts. 

Key words: galaxies: evolution - galaxies: formation - galaxies: stellar content - 
galaxies: haloes 



1 INTRODUCTION 



Ever since first introduced by IWhite fc FrenkI (119911 ). 
semi-analytic m odels of galaxy formation and evolution 



(hereafter SAM: iKauffmann et al.|ll999l: JKang et alllioosl: 
Croton et all |2006|; iBower et all l2006l: ICattaneo et all 



20061: iDe Lucia fc BlaizotI l2007l: [Monaco et all l2007l : 



Somerville et al.l l2008l : iKhochfar fc Silkl l2009l l have been 



successfully used to study how different physical processes 
determine the formation and evolution of galaxies. Based 
on halo merger trees extracted from A''-body simulations 
or analytic methods, SAMs follow the main processes 
that are thought to affect the properties of galaxies, like 
gas cooling, star formation, feedback, and merging. These 
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models provide a useful tool to study the interplay and the 
relative importance of these different physical processes. A 
de tailed review of the semi-analytic method can be found 
in lBaughl (|2006l '). 

Although various observational properties of galaxies 
are matched by the models, current SAMs still have prob- 
lems in reproducing some important observations and up to 
now, there is no semi-analytic model that is able to fit all 
the key statistical p roperties o f the o bserved galaxy popu- 
lati on. For exam ple, Guo ct ajj (|201l|) show that the model 
of iDe Lucia fc Blaizo t (i2007|) substantially overproduces the 
low mass end of the stellar mass function of galaxies. This 
problem becomes more severe as the resolution of the under- 
lying dark matter simulation increases. In addition, stellar 
mass f unctions at hig h redshifts a re normally not we ll repro- 
duced (jFontanot et al.ii200& : iMarchesini et al.ll2009l '): a good 
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match to the 2-point auto-co rrelation functio n of galaxies is 
so far difficult to achieve (e.g. lGuo et al.ll201fl ): and the rela- 
tion between the specific star formati on rate and Ralaxy stel- 
lar mass deviates fro m observations (|Somerville et al.ll2008l : 
iFontanot et al.ll2009h . 

These discrepancies may have various reasons: inaccu- 
rate physical modeling of the processes that govern galaxy 
formation; technical problems in tuning the model against a 
large set of observational constraints; the use of a fixed func- 
tional form for a poorly understood process, which overly 
limits the freedom in tuning the model; or a wrong cosmo- 
logical model adopted in the simulations. This large range 
of possibilities makes it difficult to correct identified discrep- 
ancies between m odel and ob s ervati ons. 

For example, IGuo et al.l (|201in tried to fix the over- 
prediction of t he low mass end of the ste llar mass function, 
as modeled bv lDe Lucia fc Blaizod ((20071), by increasing the 
effect of supernovae feedback. They managed to reproduce 
the amplitude of galaxy stellar mass function, and also ob- 
tained a reasonable match to the galaxy luminosity func- 
tions in different bands. However, they also predict a too 
large fraction of red galaxies at low masses, too high ampli- 
tudes of the stellar mass functions in the redshift range of 
[0.8, 2.5], and galaxy auto-correlation functions that are too 
high for galaxies less massive than 6 x I O^^Mq. M oreover, 



the high feedback efficiency as used by iGuo et al.l (120111 ) 
is physically difficult to motivate (| Benson et al.l [2003), 
and is far more efficient than various solutions adopted 
by hydrodynamical simulations ( Mac Low fc Ferraral 1 19991 : 
IStrickland fc Stevensll2000l : lAvila- Reese et al.ll201ll ). 

In this work, we therefore explore an alternative solu- 
tion. We tune the SF recipe instead of the supernovae (SN) 
feedback to study how our changes affect different statis- 
tics of galaxies, and to what degree the discrepancies men- 
tioned above can be alleviated. Most of the current SAMs 
use an analogue to the empir ical Kennicutt-Schmidt law 
(jSchmidtlll959l : iKennicuttll 19981 ) to calculate the star forma- 
tion rate (SFR) in galaxies. In this standard prescription, 
the SFR is roughly proportional to the cold gas mass, and 
scale s inversely with t he typical time-scale of a galactic disk 
(e.g. ICole et al.ll200(j ). This law is combined with a sharp 
threshold at low g as densities, below which no SF occurs 
(JKauffmannlll996l : ICroton et al]|2006l ). The use o f a such 
a threshold is motivated both th eoretically (e.g. iToomrd 
1964l:lKennicuttll989l:ISchayel2003 ') and observationally (e.g. 



Martin fc Kennicuttl[200l] ) 



This simple SF law is however likely an oversim- 
plification. In recent years observational determinations 
of the SF law in galaxies have become increasingly re- 
fined, using various gas components (HI, CO) in com- 
bination with more reliable estimates of SF, based on 
UV and IR light. These recent findings can be sum- 
marized as follows: First, there are indications that the 
threshold for SF at low mass densities is not sharp. In- 
stead, the SF effi ciency drops off as a s teep power-law at 



low gas d ensities (^Kcnnicu tt et al.l 20071: 
Wvder et al. 2009; Rovch owdhurv et al.l 



20ld ). Second, several studies find 



Bigiel et al.ll2008l: 



20091 : iBigiel et all 
that the SF rate 
is correlated more strongly with the mass of molecular 
gas ( H2) than with the the atomic gas ()Wong fc Blita 



This shows that simply correlating the star formation rate 
with the total cold gas density in models may not al- 
ways lead to realistic results. It is also not clear which 
gas mass correlates best with th e SF rate when aver - 
aging over the entire galaxy (e.g. ISaintonge et al.l [201l|). 
Third, there are both theoretical and observational indi- 
cations that the normalization of th e SF law might be 
lower at high redshift than locally (e.g. Wolfe fc Chen 20061: 
Rafelski et al.ll20ld:lG"nedin fc Kravtsovll2010l : lAgertz et aP 



201ll : lKrumholz fc DekellbOllh . This means that simply ex- 



trapolating the local relation, as done in most SAMs, may 
be incorrect. 

Se veral recent models h ave attempted to address these 



I2OO2I : iBigJel et al] I2OO8I: iLerov et al.' '200 J), probably als o 
at high redshift (|Bouche et all |2007i : iGenzel et al.ll2010l) . 



iBaugh et all (120051 ). IWeinmann et all (|2011al ). and 
iKrumholz fc Dekell (|201ll ) specifically lower the quiescent 
star formation efficiencies at high redshifts in order to bet- 
ter m atch som e obser ved properties of high redshift galaxies. 
Both lFu et ai] (|2010l ) and lLagos et all (|2010l , l201ll ) focus on 
the first two of the above points, and present SAMs with 
updated and much more detailed SF recipes in comparison 
to previous SAMs. They do not include a sharp threshold for 
SF, and their SF rate depends on the molecular gas density 
instead of on the cold gas mas s, following recent empiri- 
cal and theoretical m odels (e.g. iBlitz fc RosolowskyI 120061 : 
IKrumholz et al.ll2009l ) . Their detailed models for SF depend 
on various internal properties of the disk, like size and pres- 
sure, which are not trivial to model in a SAM. 

It is certainly worthwhile to try and include more de- 
tailed and observationally and theoretically better moti- 
vated SF law into SAMs. We present a complementary ap- 
proach in this work, without taking into account compli- 
cated processes on sub-galactic scales, like the conversion 
from atomic to molecular gas. We instead try to solve in 
a straightforward way the inverse problem, namely which 
realistic SF law at galactic scales is required to improve 
the agreement between model galaxies a nd observations. We 
start w ith a standard model, similar to lDe Lucia fc Blaizoj 
(120071 ) and Neistein fc Weinmann (2010), and assume that 
the SF rate depends on the cold gas mass in the galaxy, 
cosmic time, and in addition the host halo mass. We make 
simple change to the standard SF law that are qualitatively, 
but not quantitatively plausible, and explore how those im- 
pact on the properties of galaxies. We show that we can 
improve the agreement between SAMs and observations in 
several key aspects in this way. 

We use the model developed by Neistein & Weinmann 
(2010), a nd implement it on both the Millennium Simula- 
tion CMS. ISpringel et al.ll2005^ and the Mil lennium-II Simu- 
lation fMS-II. iBovlan-Kolchin et al.|[2009l ). The properties 
of galaxies with masses as low as ~ IO^Mq can proba- 
bly be studied r eliably with the help of these simulations 
(|Guo et al.ll201ll ). although this can be resolution depen- 
dent for extreme models (see Neistein fc Weinmann 2010). 
Our aim is to investigat e how much a model based on 
iDe Lucia fc Blaizoj (|2007l ) can be changed and improved by 
tuning the SF law alone. We thus do not change gas cool- 
ing and feedback in our models, but leave it at the default 
standard values (except in one case, for illustration, as will 
be explained). 

This paper is organized as follows. In section [2] we 
present the different models used in this work. Our start- 
ing point is a model that is based on the widely used SAM 
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of iDe Lucia fc BlaizotI ((20071) , with an improved prescrip- 
tion for hot gas stripping of satellite galaxies. We then de- 
velop four different models. Of these, models 2 and 3 include 
the key changes to the quiescent and burst mode star for- 
mation. In section [3] we present predictions for the stellar 
mass functions at low and high redshifts, the relation be- 
tween galaxy specific SF rate and stellar mass, the galaxy 
cold gas mass function, the auto-correlation functions, and 
the SF rate density as a function of redshift. A discussion 
of the implications of our results, and the conclusions are 
presented in section [l] 



2 MODELS 

The semi-analytic models presented in this paper are applied 
to both the MS and MS-II simulations. The cosmological pa- 
rameters in the simulatio ns are consistent wi th a combined 
analysis of the 2dF GRS JCoUess et al.ll200ll ') and the first 
year WMAP data (jSpergel et all 12003 ^. with f^m = 0.25, 
fib = 0.045, h = 0.73, r^A = 0.75, n = 1, and erg = 0.9. Note 
that these parameters are different from the latest WMAP 
7 year results. Both simulations follow N — 2160^ particles 
from redshift z = 127 to the present day. The MS has a par- 
ticle mass resolution of 8.6 x lO*/i~^M0, with a comoving 
box of 500 h~^Mpc on a side. The MS-II has a mass resolu- 
tion of6;£2iJ^l^lCiMai^ithabox of side 100 h'^Mpc. 

iNeistein fc WeinmannI (|2010l ') (hereafter NW2010) de- 
veloped a new formalism for modeling galaxy formation and 
evolution, which is similar to the standard SAMs, except 
that the efficiencies of processes like gas cooling, star forma- 
tion and feedback are assumed to depend only on the host 
halo mass and cosmic time. NW2010 have shown that this 
new method produces a very similar population of galax- 
ies like standard SAMs. The method is simple and flexible, 
which makes it easy to change recipes in order to fit selected 
observational constraints. 

All the models within this work are based on a sim- 
ple set of differential equations, that follow the mass of 
gas and stars with i n ga laxies. We adopt the model of 
iDe Lucia fc BlaizotI (|2007l . hereafter DLB07) as our start- 
ing point, as was don e in NW2010. The re ader is referred to 
these papers, and to ICroton et al.l (|2006l ). for more details 
on the model assumptions. Here we highlight a few features 
that will be important for the discussion below. The SF law 
is assumed to be: 



Mstar = fsiMcold - Merit) 



(1) 



fs is the SF efficiency in units of Gyr^^, and A/crit is the 
critic al mass of cold g as, below which no star formation oc- 
curs (|Kennicuttlll99a ). Satellite galaxies are followed along 
with their host subhaloes. Once the subhaloes are stripped 
and cannot be identified anymore, we compute the radial 
distance, raat, between the satellite and the central subhalo 
within the group. We then allow the satellite galaxy to spiral 
in further, and estimate the time it merges into the central 
object by using dynamical friction estimate: 



idf ~ Ctditc 



adi 



1.17KrL 



GMiln{l + Ah/Mi 



(2) 



central subhalo, Afh is its mass, and Afi is the baryonic (cold 
gas and stellar) mass of the satellite galaxy, adf describes the 
ratio of the adopted dynamical friction time over the Chan- 
drasekhar estimate. When galaxies finally merge we assume 
a SF burst of the type: 



A^star, burst — Ctburst (A-/l,cold + A-f2,cold) ; 

with 

Oburst = 0.56(Afl/Af2)°-^ . 



(3) 



(4) 



Here Afi and AI2 denote the baryonic mass in the merging 
galaxies. Following ICroton et al.l (120061 ) . this formula is de- 
rived from_^tting_the results of hydrodynamical simulations 
(e.g. ICox et al.ll2004l . see section 2.4 for more details). 

In the following subsections, we describe all the models 
used in this paper in detail. 



2.1 Model 

Our model 0) is very similar to model 0) in NW2010, and 
should thus be also similar to DLB07. For the SF law above 
(Eq. [TJ we use the same fitting functions as in NW2010: 



fs = 2.04Af?2°^nO-°-"'''"°'^'*'i^l%-'''**' , 



and 

Mcri, 



0.36 /7'Af; 



0.68,-0.52 



(5) 



(6) 



Here fc is the Chandrasekhar estimate for the dynam- 
ical friction timescale, where Vv is the virial velocity of the 



Here t is the Hubble time in units of Gyr, and A/12 ~ 
Afhaio/10^^ is the halo mass in unit of lO^^/i'^M©. 

There are a few minor modifications made here in com- 
parison to NW2010, that are related to the extended range 
in halo mass we use here (a minimum halo mass of --^ 10^ 
h"^ Mq in comparison to ~ 10^° h'^ Mq in NW2010). For 
more details on how we extend the recipes from NW2010 
to low mass haloes, the reader is referred to Appendix A. 
Fig. [T] shows that the amplitude of the low mass end of the 
stellar mass function (SMF) for model 0) is comparable to 
the DLB07 result when applied to MS-II simulation. 



2.2 Model 1 

Model 1) is the fiducial model used in this work. It is based 
on model 0) as presented in the previous subsection, but in- 
cludes two further changes: the hot gas stripping of satellite 
galaxies i s slowed down consider ably compared to DLB07 
following IWeinmann et al.l (|2010|), and a larger dynamical 
friction time is assumed. 

In the DLB07 model, the hot gas component of a galaxy 
is stripped completely once it falls into a larger group and 
becomes a satellite. Satellite galaxies subsequently consume 
their remaining cold gas due to star formation and efficient 
SN feedback, and the SF rate ceases on a short timescale 
of 1 - 2 Gyr. This leads to most satellite galaxies dis- 
playing red colour, which is in contradiction wi th obser- 
vations (|Wang et aLll2007l : IWcinmann ct al."2009J}. Follow- 
ing the model suggested bylWcimnann etal. (2 0101 ). we as- 
sume in model 1) that the hot gas component of satel- 
lite galaxies decreases at the same rate as their surround- 
ing dark matter haloes, which lose mass due to tidal strip- 
ping. This treatment provides a physically better motivated 
description of the behaviour of the hot gas component of 
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satellite galaxies and im proves agreement with observations 
(|Weinmann et al. 2010l). A simi lar model is included in the 
recent SAM of lGuo et aD (J201lh . 

The second change with respect to model 0) is that 
the parameter adf, which describes the ratio of dynamical 
friction time for galaxy mergers over the Chandrasekhar for- 
mula (Eq. [2] above), is set to 5, in contrast to the value 2 
as adopted in model 0) and in DLB07. This corresponds to 
a larger time scale for satellite galaxies to merge with the 
central galaxy, and is chosen in order to get auto-correlation 
functions in better agreement with observations for all the 
models (see below). Up to now t here is no solid consensus on 
what Odf should be in models (jBovlan-Kolchin et alj 120081 : 
IJiang et al.ll2008l : IMo et aLlbOlCl '). 

Green lines in Fig.[T]show the SMF of the fiducial model 
1), with dashed and solid lines for results when applied to 
MS and MS-II simulations. Model 1) SMF are in general 
slightly higher than model 0). This is mainly due to the 
slower stripping of hot gas in satellites adopted in model 1). 
Retaining their hot gas reservoir for longer, satellite galaxies 
can continue forming stars for a considerable time and thus 
end up with a higher stellar mass than in model 0). When 
merging into central galaxies, they also add more mass to 
their centrals. Although our choice of a larger dynamical 
friction time, with adf = 5, delays merger to some degree 
and thereby decreases the amount of mass added to centrals, 
this effect is apparently smaller than that of the modification 
to the hot gas stripping. When comparing results of model 0) 
and model 1) for the MS and MS-II simulations, the SMF 
in the MS exceed the SMF in the MS-I I at interinediat e 
masses. This is not the case for DLB07 (jGuo et al.ll2011^ . 
This excess may be related to the fact that our extension of 
the coolin g efficiencies do n ot exactly match those used for 
DLB07 in lGuo et al.l (|201ll l. 



2.3 Model 2 

In model 2), we make modifications to the SF law in the 
quiescent mode, and keep all the other components of the 
model exactly the same as in the fiducial model 1). The 
modifications include: a) Merit = 0, which means no thresh- 
old cold gas mass for SF; b) the SF efficiency depends more 
strongly on halo mass; c) the SF efficiency does not depend 
on the Hubble time. 

The SF law in model 2) is chosen such that the low 
mass end slope of SMF in the model is comparable to the 
observed SDSS result, when applied to MS simulation. This 
results in 



Merit = , 



(7) 
(8) 



Fig. [2] compares the SF efficiency (SFE, defined as the 
ratio between SF and the mass of the cold gas) in the qui- 
escent mode in model 1) and model 2), as a function of halo 
mass at two redshifts. Green lines are efficiencies in model 
1). Solid green lines are for a cold gas-halo mass ratio of 
0.04/0.09 a.t z = 0/z = 3, which corresponds to the median 
cold gas-halo mass ratio in model 1). Dashed green lines 
are for a higher cold gas-halo mass ratio of 0.13/0.14 at 
z = 0/z = 3, and dotted green lines are for a lower ratio 
of 0.007/0.05 at 2; = 0/z = 3. These values correspond to 
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Figure 1. The stellar mass functions of model 0) (black lines) 
and the fiducial model 1) (green lines). Dashed hnes are obtained 
using the MS simulation, and solid lines are for the MS-II sim- 
ulation. The red dotted line shows the stellar mass function us- 
ing the DLB07 SAM for MS-II simulation. Filled dots with error 
bars show the observed stellar mass function of SDSS DR7 result 
JLi fc Whitell2009l : IGuo et al.|[2oToh . 



the 16th and 84th percentiles of the distributions in cold gas 
masses in model 1). The SF efficiency in model 2) is shown 
as red-blue lines, which follow a power law of index 2 (plot- 
ted as black dotted lines in Fig. [2]) for haloes less massive 
than ~ lO^^'^/i~^M0. Compared to model 1), model 2) has 
a lower SF efficiency a.t z — for halo masses between 10^" 
and ~ 10"''^/i"^Mo. At z = 3, the SF efliciency is much 
lower in model 2) than in model 1) for haloes less massive 



than ~ 10^^h~^Mc. 



2.4 Model 3 

When applying model 2) to the MS simulation, the ampli- 
tude of the SMF at the low mass end decreases dramatically 
and matches the observation of SDSS (red dashed line in 
Fig. |3} . However, when applied to the MS-II simulation, the 
low mass end of SMF is still too high (red solid line in Fig.|3]). 
We have tested that even when truncating all SF in quiescent 
disk mode in haloes less massive than ~ 10 ' h~ M0, the 
problem still exists. This is because galaxy mass grows not 
only by quiescent SF in disks, but also by merger-induced 
bursts. These become more significant if the quiescent SF 
efficiency is decreased, due to the resulting increase in the 
cold gas masses. 

In model 3), we therefore further modify the SF law in 
the burst mode, while keeping the same SF law for the qui- 
escent mode as in model 2). The star burst efficiency used 
in current SAMs (Eq.|4} is derived from fitting the results of 
hydrodynamical simul ations for rnerger s of mass ratios rang- 
ing from 1 :10 to 1:1 (|Cox et al.1 12004| : iMihos fc HernqulstI 
1 19941 . [l99a) . The parameters in those simulations are set to 
make the SF in an isolated disk galaxy consistent with the 
Kennicutt-Schmidt law. As the quiescent SF of dwarf galax- 
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Figure 2. Star formation efficiencies in the quiescent mode as a function of halo mass at 2 = (left panel) and z = 3 (right panel). 
Green lines are the efficiencies in model 1). Dashed, solid and dotted green lines are for the cases when the cold gas-halo mass ratios are 
[84, 50, 16] percentile values in model 1), which are [0.13, 0.04, 0.007] at 2 = 0, and are [0.14, 0.09, 0.05] at z = 3. Red-blue line is the 
efficiency of star formation for model 2) and 3). The magenta lines are for efficiencies in model 3b). Black dotted line shows a power law 
with index 2 for reference. 



ies themselves is lower than that predicted by the Kennicutt- 
Schmidt law, it is quite possible that those simulations over- 
predict the burst efficiency for low mass galaxies. Therefore 
we introduce a halo mass dependence for the efficiency of 
merger- induced bursts in model 3), and make it inefficient 
for low mass haloes. 

For haloes less massive than Mo, we modify the burst 
efficiency to be 



Oburst = 0.56(Ml/M2)°-^x(Mhalo/Mo) . 



(9) 



while it remains unchanged for higher mass haloes. We use 
Mo — 10^^'^/i^^Mq. This critical halo mass is empirically 
determined, and roughly corresponds to the mass where the 
relation betwe en galaxy stellar mass and halo mass ch anges 
its slope (e.g. IWang et all l2006l : iMoster et all I2OI0D. and 
galax y formation efficiency reaches its maximum ( Guo et al.l 
[2OI0I). Physically, SN feedback and reionization are believed 
to cause a low galaxy formation efficiency at low halo masses, 
and AGN feedback may be responsibl e for the low efficiency 
in high mass haloes (jGuo et al.ll2010l ). The combined effect 
of these mechanisms may be weakest at this critical halo 
mass, explaining the peak in galaxy formation efficiency in 
the current theory. Therefore we choose this value as the 
threshold below which the SF efficiency is assumed to be 
suppressed. The typical stellar mass of galaxies that reside 
in those haloes is about 10^"' ^Mq, below which model 0) 
and model 1) predict too many galaxies and thus too much 
SF. By decreasing the SF for galaxies with halo mass less 
massive than this critical value, the amplitude of the low 
mass end of the stellar mass function can be suppressed ef- 
fectively. The functional form in Eq. [9] is chosen in order 
to fit the amplitude of the SMF at the low mass end when 
applied to both MS and MS-II simulation. 

We have also tested a model in which we shut off all SF 
in the burst mode in low mass haloes and keep the SF in the 
quiescent mode the same as in the fiducial model 1) . For such 
a model, the low mass end of the SMF is still higher than ob- 



servation. This indicates that SF in both the quiescent and 
the burst mode must be modified simultaneously to suppress 
the numbers of low mass galaxies effectively. Note that the 
values of the power law indices that determine the depen- 
dence of SF on halo mass for quiescent and burst modes 
might have some degeneracy. However, we do not study the 
degeneracy of the two modes of SF in this work any further, 
but focus on the qualitative effect of modifying SF in each 
mode. 



2.5 Model 3b 

In model 3b), we test the effect of including a dependence 
of the SF efficiency on Hubble time. Model 3b) is almost 
identical to model 3), except that the SF efficiency in the 
quiescent mode is assumed to depend on cosmic time, in the 
same way as for the fiducial model 1). The SF efficiency is 
rescaled to match again the low mass end of SMF in SDSS 
observation, resulting in: 



/. 



1.74MiV*10 



-0.30[logJ\/i2l^ ,-0.82 



(10) 



The SF efficiency in model 3b) is also shown in Fig. [2] It 
is higher than the SF efficiency in model 3) at high redshift, 
and is lower at 2; = 0. We will see in the next section that 
compared with model 3), model 3b) results in a similar SMF 
at z = 0, and higher amplitude of SMF at higher redshifts 
of z '^ 2 — 3. The lower SFR at redshift results in low 
mass galaxies with a lower SSFR than in model 3). Besides, 
model 3b) predicts a higher SFR density at redshifts higher 
than ~ 2, and a lower SFR density at lower redshifts. The 
general results are comparable with model 3) though. 

2.6 Model 4 

With a modified SF law in low mass haloes, model 3) and 
3b) are able to reproduce many observational statistics of 
low mass galaxies, as we will show in in Sec.|3] However, the 
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properties of high mass galaxies differ significantly from ob- 
servations. The most obvious deviation is that the modeled 
massive galaxies are in general too active. This is mainly 
due to the inclusion of slower stripping of hot gas for satel- 
lite galaxies in our models. The deviation can be alleviated 
by allowing for less efficient cooling in massive haloes. Phys- 
ically, this corresponds to mechanisms like stronger AGN 
feedback effect that prevents gas from cooling. 

In model 4), we apply further modifications of cooling 
and SF to model 3), focusing on massive galaxies. Model 4) is 
presented as a simple test to see if the properties of massive 
galaxies can be better fitted, while keeping the the treatment 
of low mass galaxies unchanged. The modifications include: 

• Lower cooling efficiencies are assumed for haloes more 
massive than 1O^^'^'''/!.^^M0, as shown in Fie:. IA2l in the Ap- 
pendix. 

• SF in both the quiescent and the burst mode is stopped 
completely in haloes more massive than 5 x lO^^/i^^Me at 
z< 1.3. 

• The dynamical friction time is assumed to be depen- 
dent on Hubble time and is shorter at higher redshift, with 
adf = 5 X (t/13.6)°'^ i nstead of an = 5 in model 3) (see 
IWeinmann et al.|[2oTlal ). 

These modifications are done to fit the observed prop- 
erties of massive galaxies better. The first two modifications 
make massive galaxies much more passive than in model 3). 
T he change in the dynami cal friction time follows the idea 
of IWeinmann et al.l (|2011al ). As discussed there, the merger 
time in the standard model may be overestimated by an or- 
der of magnitude at high redshift, mainly due to the more 
radial orbits of high redsh ift satellite galaxies (jDekel et al.l 
I2OO9I : iHopkins et al.l I2OI0I ). Therefore we assume a time- 
dependent dynamical friction time, which gives a better fit 
to the SFR density and SMF at z > 2, and does not affect 
the other statistics studied in this work much. 



3 RESULTS 

In this section, we show statistical results for the galaxy 
population produced by the models presented in the last sec- 
tion, including galaxy SMF at both low and high redshifts, 
the specific SF rate-stellar mass relation, the cold gas mass 
function, the projected two point correlation functions, and 
SF rate density as a function of redshift. We compare these 
results with observations, and analyze the effect of different 
modifications of SF laws on those statistics. 



3.1 Stellar mass function at 2 = 

The SMF at 2 = is the main quantity that was used to 
constrain the parameter values in our models, a nd is plot- 
ted in Fig. [3 Similar to previous S AMs (DLB07, IGuo et all 
I2OIII '). the fiducial model 1) predicts too many low mass 
galaxies. As a result of the modifications in the SF law in 
low mass haloes, model 2) applied to the MS fits the SMF 
for stellar masses less massive than lO^'^M©. However, when 
applied to the MS-II, model 2) still predicts too many low 
mass galaxies. By suppressing the SF in the burst mode in 
low mass haloes, model 3) gives a reasonable fit to the ob- 
served SMF, for both MS and MS-II simulations. Including 
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Figure 3. Stellar mass functions at 2 = 0. Green, red, blue, 
magenta and cyan lines are for model 1), 2), 3), 3b) and 4) re- 
spectively. Solid and dashed lines are for MS-II and MS results. 
Black symbols with error bars are based on SDSS DR7 includ- 
ing both counting and cosmic variance uncertainties I ILi fc White! 
|2009I : IGuo et al.ll2010l). The grey diamonds show the functional fit 
to the SMF from lBernardi et al.l 1I2OIOI) , for galaxies more massive 
than IQl^'^Mo. 



a dependence of the SF on cosmic time in model 3b) does 
not affect the SMF noticebly. For model 4) which changes 
further the cooling and SF in high mass haloes, the SMF 
at z=0 is somewhat lower at the massive end, and is more 
consistent with observations. 

As mentioned in sec. (2] we have tested that when shut- 
ting off all SF in either quiescent mode or burst mode in 
haloes less massive than 10^^'^/i~^Mq, the predicted SMF 
still exceeds observations at low mass end, since SF from 
the other mode compensates. This indicates that SF in both 
quiescent and burst mode must be modified simultaneously 
to fit the observed SM F, as is don e in m odel 3) and 3b). 
This also explains why iLagos et al.l (|2010l ) find no obvious 
change in the resulting SMF when only modifying the SF in 
the quiescent mode in their SAM. 

3.2 SSFR-Mstar relation 

The specific SF rate (SSFR) is defined as the ratio between 
the SF rate and the stellar mass of a galaxy. The relation 
between SSFR and galaxy stellar mass is a fundamental 
observable which needs to be reproduced by a successful 
model. Observationally, there is a clear trend that high mass 
galaxies are passively evolving with low SSFR, while low 
mass galaxies have high values of SSFR (|Salim et al.l 120071 : 
ISchiminovich et al.l 120071 ). SAMs, however, usually predict 
similar, if not lower SSFR in low ma ss galaxies than that in 
massive ones (jFontanot et al.l 120091 1. Comparisons that fo- 
cus on galaxy colours also show this discrepancy, with low 
mass galaxies having redder colours than obs erved (Quo et 
al. 2011, but see also IWeinmann et al.|[2011bl ). 

Fig. |4] shows the SSFR-Mstar relation in our models, 
compared with SDSS observational result. The upper panels 
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Figure 4. Contour plot of SSFR— Mstar relation in different models compared with observation. Upper panels are for models applied to 
the MS and lower panels are for MS-II. Blue dashed contours are SDSS DR7 observation, derive d from the MPA-JHU re lease of the SPSS 
DR7 catalogue, with stellar m asses based on fits to the p hotometry following the philosophy of lKauffmann et al.l 1120031 ) and lSalini et al.l 
112007]), and SF rates based on iBrinchmann et al.l l|2004l ). Black contours are the results of models 1), 2), 3), 3b) and 4). When plotting 
the model results, we set the minimum SSFR value to be 10~*Gyr~^. Contours levels show regions that enclose 38, 68, 87 and 95 percent 
of galaxies. 



are for models applied to the MS simulation, and the lower 
panels show the MS-II results. In each panel, black con- 
tours are the model results, while blue contours are SDSS 
results. Observationally, galaxies reside in two distinct se- 
quences in the SSFR-stellar mass plot: an active sequence 
with high SSFR that is more prominent for low mass galax- 
ies, and a passive sequence with lower SSFR that contains 
mainly massive galaxies. The observational location of the 
passive sequence, however, is not well determined, due to 
the uncer tainty on measur ing the SSFR for galaxies with 
little SF (|SaUm et al.ll2007l ). In addition. Fig. [5] shows the 
median SSFR as a function of galaxy stellar mass. Colored 
lines are predictions from our models combined with MS-II. 
The black solid line is the result of SDSS DR7 galaxy sam- 
ple, which shows clearly that more massive galaxies have in 
general lower SSFR. 

Comparing the model results applied to the MS and 
MS-II simulations respectively in Fig. 4, it is obvious that 
the resolution of simulation has a large effect on the mod- 
eled SSFR for low mass galaxies, especially for models 1) 
and 2). This highlights the fact that resolution can have a 
large impact on semi-analytical model predictions (see also 
Guo et al. 2011). We note that the resolution dependency is 
weaker for models 3)-4), only affecting the passive sequence 
of galaxies whose location is anyway uncertain (see above). 
To compare the model prediction with observation for low 
mass galaxies, we focus on results of the models combined 
with the MS-II simulation. Similar to previous SAMs, model 
1) predicts more passive low mass galaxies than observed, 
and the active sequence in model 1) is lower than observed. 
Its predicted median SSFR is lower than the SDSS result for 
galaxies less massive than lO^'^M©. When modifying the SF 
in the quiescent mode in low mass haloes in model 2), more 
low mass galaxies become passive, and the median SSFR 
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Figure 5. The median SSFR as a function of galaxy stellar mass. 
The black solid line shows observational results from the SDSS 
DR7 galaxy sample. Green, red, blue, magenta, and cyan lines 
are predictions from our models 1), 2), 3), 3b) and 4), combined 
with MS-II. 



is even lower than in model 1). This is due to the in gen- 
eral lower SF efficiency assumed in model 2) for haloes less 
massive than ~ 10^^'^h~^MQ, as shown in Fig. [2] 

Suppressing the SF in the burst mode in low mass 
haloes, model 3) and 3b) predict a lower fraction of pas- 
sive galaxies, with higher median SSFR for low mass galax- 
ies. This is because the lower burst efficiency leaves galaxies 
with more cold gas by 2 = 0. The results are then in reason- 
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Figure 6. Cold gas mass functions at z=0. Green, red, blue, 
magenta and cyan lines are for model 1), 2), 3), 3b) and 4), for 
MS-II results. Cross symbols show the Schechter function fit for 
the co ld gas mass function according to lObreschkow &: Rawlingg 
||2009| '). 



able agreement with observations. When comparing these 
two models in detail, model 3b) predicts on average more 
passive low mass galaxies and is therefore a less good match 
to observations. The location of the active sequence is lower 
than in observation, and thus not as well reproduced as in 
model 3). 

For all models l)-3b), the SSFR for massive galaxies are 
significantly overpredicted by the model. This is corrected in 
model 4), by modifying the cooling and SF in massive haloes. 
Overall, model 4) is in good agreement with the observed 
SSFR-Mstar relation at all stellar masses. 



3.3 Cold gas mass function 

Fig. [5] presents the cold gas mass functions at z=0. Green, 
red, blue, magenta and cyan lines are for model 1), 2), 
3), 3b) and 4), using the MS-II results. Crosses show the 
Schechter function fit to the cold gas (HI+Hg+heliu m) mass 
function given in lObreschkow fc Rawlingsl (|2009|), which 
is deriv ed from the HI m ass function of HIPASS obser- 
vation (jZwaan et al.l [2005|), combined with a theoretical 
m odel to describe the H2/ HI mass ratio. As pointed out 
in lObreschkow et al.l ((20091), semi-analytic models only dis- 
tinguish between the hot and the cold gas phase. How- 
ever, observations have revealed the existence of a signif- 
icant amount o f warm and ioniz ed gas in the Milky Way 
(JRevnoldsll2"00i l. This motivated fobreschkow et all (J2009i l 
to devide their model cold gas mass function by a factor 
of 1.45. We follow their example here, but we stress that 
this factor is obviously uncertain and may well be higher, 
depending on galaxy mass and redshift. 

The cold gas mass function of the fiducial model 1) ex- 
ceeds observations. As e xpected, it is somewhat h i gher than 
for t he DLB07 model (jObreschkow et al.1 l2009l : IFu erall 
I2OI0I ) , because of the slower stripping of hot gas for satellites. 



The other models predict even higher amplitudes. With on 
average lower SF efficiencies in low mass haloes in model 2), 
3) and 3b), there is much more cold gas left in the galax- 
ies. This is not surprising, since we have left cooling and SN 
feedback efficiencies at the default DLB07 values. Model 4) 
predicts similar result of cold gas mass function as in model 
3). 

Although there are some uncertainties related to the 
determination of cold gas mass function in observation 
dKeres et al.ll2003l : [Obreschkow fc Rawlingsl I2OO9I : IFu etlo] 
l2010f ). it seems that our models with modified SF predict 
drastically too much cold gas mass in galaxies. Thus, while 
the modifications to the SF law we suggest here improve the 
agreement between the statistical properties of the galaxy 
population in many aspects, our simple model is clearly not 
the final answer to all the problems that exist. This will be 
discussed in more detail in section 4. 



3.4 Correlation function 

The two point auto-correlation function of galaxies is a fun- 
damental measure of the spatial distribution of a certain 
population of galaxies an d is well determined by observa- 
tions fe.g. iLi et al.ll2006l ). It is howeve r not w ell fitted by 
current SAMs. For example, IGuo et al.l (|201ll ) over-predict 
the correlation functions for low mass galaxies and on small 
scales. They suggest that the over-prediction of clustering of 
galaxies on small scales in the models can be explained by 
the too high value of as adopted by the Millennium sim- 
ulation, which is 0.9 compared to . 81 su ggested by the 
WMAP-7 year resuh (JKomatsu et aLlbOllI ). However, it is 
not certain that the Guo et al. SAM combined with the 
correct cosmology would gi ve correlation func tions in agree- 
ment with observations (see I Wang et al.ll2008l ). Here we test 
the effect of modifying the SF law on the resulting corre- 
lation f unctions us i ng the same underlying A'^-body simula- 
tions as IGuo et aD (|201ll l. 

In our models, galaxy positions are determined by the 
positions of haloes/subhaloes they reside in. For galaxies 
that have lost their host subhalo due to stripping, we use 
the location of the most-bound-particle of the last identified 
subhalo. Since we use the same dynamical friction prefactor 
of Qdf = 5 for models 1), 2), 3), 3b) and 4), galaxy locations 
in all these models are exactly the same. The only differences 
are in the stellar masses of galaxies that vary due to the 
different SF laws used. Also, the number of galaxies might 
be slightly different due to the effect of the stellar mass on 
the dynamical friction time. Consequently, the differences in 
the correlation functions between the models are mainly due 
to the different stellar mass assigned to each galaxy. 

Fig. [7] shows the projected correlation functions of 
galaxies in di fferent stellar mass bins, computed in the 
same way as in lNeistein et al.l (|201ll ). Model 1) over-predicts 
the correlation functions for galaxies less massive than 
10"'^^Mq, and at scales smaller than ~ l/i"^Mpc. With 
a modified SF law in low mass haloes, correlation functions 
become lower at small scales for low mass galaxies, which 
brings model 3) into agreement with the observational re- 
sults at all stellar masses. The success of model 3) in repro- 
ducing the correlation functions shows clearly that changes 
in the SF law, that lead to changes in the relation between 
stellar mass and subhalo mass, can significantly affect cor- 
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Figure 7. Projected two point auto-correlation functions of galaxies in diffe rent stellar mas s bins. Circles with error bars are SDSS DR7 
results ||Guo et al.ll201M . calculated with the same method as presented in iLi et al.l l|2006h . Green, red, blue, magenta, and cyan lines 
are for model 1), 2), 3), 3b), and 4) respectively. Dotted (dashed) lines are for models based on the MS-II (MS) results. For stellar mass 
bins with log(Mstar) > 9.77, only MS results are shown. For the lowest mass bin, only MS-II results are presented. 



relation functions. Thus, it is in general very hard to say 
whether a mismatch with observed correlation functions in- 
dicates a problem with baryoni c recipes, or with c osmology. 
We thus confirm the results by I Wang et al.l (|2008l 'l who find 
that a similar match to observed correlation functions can be 
obtained by SAMs using different cosmological parameters, 
depending on the detailed baryonic recipes. 

Model 3b) gives similar results to model 3) for galaxies 
more massive than IO'^'^'^^Mq. For lower mass galaxies, the 
correlation functions of model 3b) are higher than for model 

3) on small scales. Correlation functions of galaxies in model 

4) are in general similar to model 3), with a somewhat worse 
fit in the stellar mass bin of logMstar = [10.77, 11.27]. 



3.5 Stellar mass functions at high z 

At redshifts less than ~ 0.8j_ current SAMs like 



[Monaco et"all (|2007l ') and ISomerville et all (|2008l ') predict 
SMF consistent with observations. At higher redshifts, 
however, models normally over-pred ict the abundance of 
galaxies less ma ssive than lO'^'^Mo (jFontanot et al.l |2009| : 
IGuo et al.l I2OIII ) . The same is true fo r the K-band lumi- 
nosity function (jHenriaues et aLluOllI ). This may indicate 
that SF at redshifts above 0.8 is not modeled correctly in 
these models. 

Fig. [S] gives the results of SMF at higher redshifts in 
different models studied in this work. Model results are pre- 
sented at redshifts of 0.8, 1.2, 2.0, 3.0, and are convolved 



with a Gaussian error of deviation 0.25 dex in logMstar, to 
account for t he various errors in es timating stellar masses in 
observations (|Fontanot et al.ll2009l ). Observation results are 
shown at comparable redshift ranges to the models. Gold 
symbols indicate data points below the limiting stellar mass 
of the observed galaxy samples, whe re incompleteness could 
be significant (|Pozzetti et al.l [20071 : Ka iisawa ct al. 200l|). 
All observational s tellar masses of galaxies are normalized to 
the Chabrier IMF (|Chabrieijl2003l '). to be consistent with the 
previously shown SDSS SMF, and the model derived SMF. 

The fiducial model 1) predicts too many low mass galax- 
ies at all redshifts. With a modified SF law in both modes, 
model 3) and 3b) both give consistent SMF with observa- 
tions up to redshifts of around 3. Differences of these two 
models can be seen at redshifts higher than 2, where model 
3b) predicts a SMF with a higher amplitude. This is because 
the SF efficiency in model 3b) depends on time and is higher 
at high redshifts. 

In model 4), with the assumption that the dynamical 
friction time is shorter at high redshifts, SMF at redshifts 
of around 2 — 3 are higher than in model 3), and some- 
what closer to observations. The SMF in model 4) happen 
to be quite similar to the results of model 3b). This refiects 
the degeneracies inherent in SAMs, in this case between the 
dependence of dynamical friction time on redshift, and the 
dependence of SF rate on redshift. 
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Figure 8. Stellar mass functions at rcdshifts above zero. Green, red, blue, magenta and cyan lines are for model 1), 2), 3), 3b) and 4). 
Mode l results combined with MS-II are shown at redshifts of 0.8, 1.2, 2.0, 3.0. Black points are obs e rvatio nal results from lKaiisawa et al.l 
1120091 ). for all four redshift bins. In the upper left panel, observational results from lPozzetti et al.l l|2003) at 0.7 < z < 0.9 are plotted as 
diamonds. Gold s ymbols are the da t a poin ts below the completeness limit in these two observational samples. In the lower two panels, 
observations from iMarchesini et al.l 1 120091) fo r 1.3 < 2 < 2. and 2.0 < z < 3.0 are shown as diamonds. Observed stellar masses of 
galaxies are normalized to the Chabrier IMF llChabrieill2003l) . 



3.6 SFR density 

Fig. |5] shows the star formation rate (SFR) density as a 
function of redshift, as predicted by different models com- 
bined with MS si mulation. Black crosses are observational 
results compiled bv lHopkinsI (|2007l ). The grey shaded region 
sh ows the 1-a confiden ce level of the observational results 
by IWilkins et al.l (|2008h , which are derived indirectly from 
the evolution of the stellar mass fu nction. Gold symbols are 
the results of lBouwens et al.l (|2009l ). including the contribu- 
tions from highly dust obscured galaxies and ULI RGs. Stel- 
lar m asses are normalized to the Chabrier IMF (jChabrieil 
I2OO3I ). 

The SFR density of the fiducial model 1) shows a contin- 
uous increase with redshift and peaks at a redshift of around 
3 — 4, which is clearly a higher redshift than in observations. 
This offset is similar to the one present in the SAM of Guo 
et al. (2011). With a modified SF efficiency in both the qui- 
escent and the burst mode, the SFR density of model 3) and 
3b) drops dramatically at higher redshifts than z ~ 2. Due 
to the dependence of SF on Hubble time assumed in model 
3b), this model gives a higher SFR density at high z, and a 



lower SFR density at low z than model 3). At low redshifts, 
both model 3) and 3b) predict SFR density higher than the 
fiducial model 1), which lies slightly above the observational 
values. 

With further suppression of cooling and SF in massive 
haloes in model 4) , the observed sharp decline of SFR den- 
sity towards low redshifts appears. Predictions of model 4) 
are within the observational constraints, whil e the SFR den- 
sity p eaks at redshift of around 2. Recently. iMagnelli et al.l 
(|201lh studied the evolution of dusty infrared luminosities 
function using Spitzer data. Assuming a constant conversion 
between the IR luminosity and SFR, they find that the SFR 
density of the Universe strongly increases towards z — 1.3, 
and stays constant out to z = 2.3. Model 4) matches the 
result of their observation. 



4 DISCUSSIONS AND CONCLUSIONS 



We u se the method developed by iNeistein fc WeinmannI 
(J2OIOI), combined with both the Millennium (MS) and 
Millennium-II (MS-II) cosmological simulations, to study 
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Figure 9. Cosmic star formation density as a function of rodshift. 
Green, red, blue, magenta and cyan lines are for models 1), 2), 3), 
3b) and 4) respectively, based on the M S simulation. Bl ack crosses 
are observational estimates compiled bv lHopkinsI l|2007h . The grey 
shaded region shows th e \-a confidenc e level of the observational 
result, as co mpiled by [ Wilkins et al.l 1(2009). Gold symbols arc 
the results of iBouwens et al.l 1 120091 ). including the contributions 
from highly dust obscured galaxies a nd ULIRGs. Stellar masses 
are normalized to the Chabrier IMF jChabricr., 20031 ) . 



the effect of modffying the star formation (SF) recipe in 
low mass galaxies. We show that by modifying SF in both 
the quiescent and the burst mode, the stellar mass function 
observed in the local Universe can be reproduced well down 
to IO^'^Mq. Simultaneously, the models can fit the observed 
median SSFR-Matar relation for galaxies less massive than 
IO^'^Mq, the correlation functions for galaxies more massive 
than IO^'^^Mq, the stellar mass functions up to redshift of 
around 3, and the general trend of the SFR density as a 
function of redshift. 

The modifications to the SF recipe in our models with 
respect to standard SAMs (e.g. DLB07) include: 

(i) no sharp threshold in the cold gas mass for SF; 

(ii) letting the SF efficiency in the quiescent mode depend 
on host halo mass; 

(iii) removing the dependence of the star formation rate 
on Hubble time (which came via the dependence on disk 
dynamical time); 

(iv) a lower star burst efficiency in low mass haloes; 

(v) additional modifications of cooling and SF in massive 
haloes to match the properties of high mass galaxies. 

Model 2) only includes the changes to the quiescent 
mode of SF, i) - iii). We show that this is not enough to re- 
produce the low mass end of the SMF, as star formation in 
the burst mode compensates for the changes to t he SF law in 
the qu iescent mode. This is also the reason whv lLagos et al.l 
((20101) have found that their changes to the SF law in the 
quiescent mode does not change the resulting SMF much. 
In model 3), we have thus additionally decreased SF in the 
burst model (modification iv), which results in a clearly im- 
proved SMF. We also considered a model 3b) which is sim- 
ilar to model 3), except that we allow for the usual time- 
dependence of the SF law (i.e. do not make modification 



iii). Results of model 3b) are similar to model 3) out to 
2 ~ 2, except that the SSFR-Mgtar relation is slightly less 
well reproduced. It is thus not clear whether modification 
iii) is necessary. 

Note that removing the time-dependence of the star 

formation efficiency is a significant change with r espect to 

revious models. Even i n the recen t models of [Fu et al.[ 



20101 ) and JLagos et al] (|2010l . 1201 ih . where SF efficiency 
does not depend on the disk dynamical time of galaxies, a 
time-dependence enters via the conversion efficiency from 
atomic to molecular gas, that depends on the gas density. 
In order to justify the behaviour we suggest in model 3), we 
would need to postulate a mechanism that scales with time 
in the opposite way than usually assumed, like for example a 
mctallicity-dcpcndcnt conversion of atomic to molecular gas 
(Krumholz & Dckcl 2011,). We note that a weak dependence 
of the SF efficien cy on cosmic time was seen in hydrodynam- 
ical simulations ((Neistein et al.|[2011bl ). although the reason 
for this behaviour is still not clear. 

In our models 2) - 4), the ratio between the quies- 
cent star formation rate and the cold gas mass, M^tar / Mcoid 
which is equal to the gas consumption timescale, is roughly 
proportional to M^^i^ for low mass haloes and is almost in- 
dependent of halo mass for high masses. This m ay in fact 
be su pported by recent observations. Recently, [Shi et al.[ 
(|201lh derive an extended Schimdt law from an observed 
galaxy sample that extends over 5 orders of magnitude in 
stellar density, including galaxies with low surface bright- 
ness. They find that Afstar/A/coM is proportional to Mftf^^, 
with a 1-cr scatter of 0.4 dex. The stellar mass of galax- 
ies has been claimed to obey a tight r e lation with the host 
halo mass ((Conrov fc WechsleJ |2009| : JMoster et all l2010l : 
IGuo et al.ll2010l ). It is proportional to M^„f„ for low mas s 
haloes and to M^'^^ at high mass end (jWang et al.ll2006l ). 
Without considering the scatter of the relation, this indi- 
cates that the observed Matar/Afcoid is roughly proportional 
to Mjjg^o at low masses and M^^J^ at high mass end. The 
dependence on halo mass in our models is therefo re quite 
close to the observational result bv lShi et al.l (|201l( ). 

In all our models, the cold gas mass function of galax- 
ies is dramatically over-predicted. This is because in DLB07, 
the total amount of cold g as and stellar mass exceeds the ob- 
serve d total amount (see lObreschkow et al.l l2009l : IFu et al.l 
I2OI0I ). Decreasing SF rates in low mass galaxies while let- 
ting cooling and feedback recipes remain unchanged natu- 
rally results in an overproduction of the cold gas mass in 
low mass galaxies. This is not only a prob lem for the mod - 
els we present here. As shown recently bv lLu et al.l (|201H ). 
when the model K-band luminosity function is forced to 
fit the data in the local Universe, the cold gas mass func- 
tion is dramatically over-predicted in all the semi-analytic 
models they study. This is consistent with our models over- 
predicting the cold gas mass functions when we fit the stel- 
lar mass function at z=0. On the other hand, with similar 
cooling and supernovae feedback recipes, the models that 
do fit the cold gas mass function in turn h ave p roblems in 
reproducing the stellar mass functions (e.g. IFu et a l.„ 2010l; 
lObreschkow et al.l 120091 ). In the meantime, approaches like 
increasing the SN feedback, that suppress the total amount 
of cold ga s and stellar ma ss, lead to several other serious 
problems (JGuo et al.l 120111 ). These results show again that 
it is currently difficult for a single model to fit all observa- 
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tions, as mentioned in Sec.[Tl unless we allow for free tuning 
of all recipes, including cooling (see NW2010). 

Although the modifications of the SF law presented in 
this work help to improve the agreement with several ob- 
served statistical properties of galaxies, and seem to follow 
a similar scaling like the observed SF efficiencies in galax- 
ies, the normalization of the SF efficiency cannot be correct. 
This is clear from the fact that our models overpredict the 
cold gas mass function. With a lower cold gas fraction in 
galaxies, the SF efficiency would obviously have to be higher 
than currently assumed in the models, to obtain the same 
SF rate. 

In Fig. IIOI we compare the SFR-HI mass relation in 
model 3) combined with MS-II (right panel), with a re- 
cent observation of the FUV derived SFR-HI ma ss relation 
for g alaxies within ~ 11 Mpc of the Milky Way (JLee et al.l 
I2OIII . left panel). The cold gas mass in model 3) is con- 
verted to the HI mass to be compared with observation, us- 
ing a correction including three factors. First, as for Fig. |6l 
the cold gas mass from model 3) is divided by a factor of 
1.45, to account f o r a w arm ionized gas phase, as done in 
lObreschkow et al.l ()2009l ). Second, we multiply the results 
by a factor of 0.76 to remove the con tribution of helium and 
heavier elements (jPower et al.ll201Cl ). Fi nally, the hydrogen 
gas mass is divided by 1.4, as adopted bv I Power et al.l (|2010l ) 
and lLu et al.l (|201 j ) , to remove the contribution of H2 ■ Since 
the galaxies in the observed sample have stellar masses less 
than ~ IO^^Mq, we present model results for galaxies with 
8 < log(Mstar) < 10. The median relation from the observa- 
tions is pl otted as re d line i n each panel. Although the obser- 
vations of I Lee et al.l (|201H ) are limited to a small volume of 
space, the obviously different relations in observations and 
in model 3) indicate that the SF efficiency in model 3) is 
indeed much lower than in reality. 

Fitting the most important observed properties of 
galaxies is not a trivial task. It is not clear up to now if 
this difficulty in modeling galaxy properties reflects a fun- 
damental problem in our understanding of the dark matter 
universe, or if it is mainly due to an insufficient understand- 
ing of the baryonic physics involved in galaxy formation and 
evolution. For example, perhaps SAMs miss an important in- 
gredient of galaxy formation, like a mechanism that preheats 
the gas in the universe so that it c annot cool to low mass 
haloes (|Mo et al.ll2005l . but see also ICrain et al.l (|2007l ')'). or 
a form of feedback t hat mainly heats low entropy gas in 
high redshift haloes (JMcCarthv et al.|[201ll '). Alternatively, 
it could be that cooling is over-efficient in the current SAMs 
for some reason. 

With the tests carried out in this work, we show that 
only modifying star formation in the SAMs can already 
improve agreement with observations in several important 
aspects. Up to now, a high-resolution SAM that matches 
the SSFR-stellar mass relation at z = 0, the SMF and 
its evolution, the correlation functions and the cold gas 
fractions of galaxies simultaneously, does not yet seem to 
exist. To find such a model, approaches tha t allow scan- 
ning of a large parameter spac e (|HcnriQu cs et al.l |2009| : 
ILu et al.ll20ld : iBower et al.ll2010l ). and approaches that al- 
low deviations from the usually assumed functio nal forms 
for physical recipes (jNeistein fc Weinmannll2010l ). or that 
include other physical proce sses than currently considered 
(JHenrigues fc ThomasI I2OI0I ) may be promising. However, 



even if one or several models are found that do indeed repro- 
duce all these fundamental observables, it will be important 
to identify degeneracies, and to verify whether the models 
are physically plausible and can be brought into agreement 
with alternative approaches, like predictions from hydrody- 
namical simulations. 
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APPENDIX A: COOLING EFFICIENCIES 

Our model 0) is similar to model of NW2010, with a few 
modifications to the cooling efficiencies, which we explain 
below. Model of NW2010 is adapted to the MS simula- 
tion, and thus only includes efficiencies down to the resolu- 
tion limit of this simulation, which is ~ W^'^h'^Mir, in halo 
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Figure Al. Cooling efficiencies as a function of halo mass and 
cosmic time in model 0), the fiducial model 1), and also in model 
2) and 3), 3b). The gray scale shows Log values of the cooling 
efficiencies in units of log[Gyr~^]. X-axis is the cosmic time, with 
13.7 Gyr corresponding to the present day. 



mass. To apply it to the MS-II simulation, we need efficien- 
cies down to a lower halo mass of ~ W^h~^h~^MQ. This is 
straightforward for most processes, as they are parametrized 
by functional forms which can easily extended to lower halo 
masses. The only exception are the cooling efficiencies, /c, 
defined as AtticooI ~ /ci^hotAi, where Anicooi is the amount 
of gas that is cooled within a time-step At, and TTihot is the 
mass of hot gas. 

In DLB 07, the treatment of gas cooling follows the de- 
scription of ICroton et al.l (120061'). where coo l ing effi ciencies 
are calculated according to IWhite fc FrenkI (|l99ll ) , assum- 
ing an isothermal gas density profile. In model of NW2010, 
cooling efficiencies are median values computed from a large 
statistical sample of galaxies in DLB07, for each bin of halo 
mass and cosmic time. The tabulated values as a function 
of halo mass and time in NW2010 follow no specific func- 
tional form and can thus not easily be extended to lower 
mass haloes. The obvious solution would be to extract those 
values from the DLB07 SAM as applied to MS-II, but this 
is not possible for technical reasons. Therefore, we estimate 
cooling efficiencies for low mass haloes such that (i) the gen- 
eral trends of cooling efficiency as a function of halo mass 
and redshift are preserved, and (ii) the resulting low mass 
end slope and amplitude of stellar mass function at z=0 is 
similar to the result of DLB07 when applied to MS-II Siin- 
ulation (see Fig.[T)). In this way our estimates of the cooling 
efficiencies at a given halo mass should be similar as in the 
DLB07 model, when averaged over all redshifts. However, 
we note that the cooling efficiencies in a given redshift and 
halo mass bin may differ. Apart from extending the cooling 
efficiencies to lower mass haloes, we also apply some smooth- 
ing to the original cooling efficiencies found by NW2010, in 
order to smoothen the stellar mass functions. 

Fig. lAll shows cooling efficiencies as a function of halo 
mass and cosmic time in model 0). These values are given 



Table Al. Values of cooling efficiencies in units of Gyr~^. The 
values shown here are identical for models 0), 1), 2), 3) and 3b), 
and arc plotted in Fig. lAll Halo mass is in units of h~^yiQ and 
is shown at the left column, time is in Gyr. 



Log(Mh„;o) 


t = 0.80 


2.24 


3.38 


5.97 


10.27 


13.58 




z = 7 


3.0 


2.0 


1.0 


0.3 





8.00 


-1.30 


-1.60 


-2.00 


-2.00 


-2.00 


-2.00 


8.25 


-0.90 


-1.50 


-1.90 


-2.00 


-2.00 


-2.00 


8.50 


-0.80 


-1.30 


-1.50 


-1.60 


-1.80 


-2.00 


8.75 


-0.50 


-1.00 


-1.30 


-1.50 


-1.50 


-1.80 


9.00 


-0.30 


-0.80 


-1.10 


-1.30 


-1.50 


-1.60 


9.25 


0.10 


-0.80 


-1.10 


-1.30 


-1.30 


-1.50 


9.50 


0.80 


-0.50 


-0.80 


-1.10 


-1.20 


-1.30 


9.75 


1.30 


-0.10 


-0.50 


-1.00 


-1.10 


-1.20 


10.00 


1.30 


0.20 


-0.30 


-0.70 


-1.00 


-1.10 


10.25 


1.53 


0.50 


-0.07 


-0.40 


-0.76 


-0.93 


10.50 


1.36 


0.81 


0.29 


-0.23 


-0.57 


-0.76 


10.75 


1.32 


0.83 


0.69 


-0.06 


-0.38 


-0.53 


11.00 


1.13 


0.42 


0.42 


0.41 


-0.32 


-0.49 


11.25 


0.77 


0.33 


0.17 


-0.06 


-0.28 


-0.44 


11.50 


0.71 


0.21 


0.02 


-0.32 


-0.41 


-0.67 


11.75 


0.43 


-0.15 


-0.45 


-0.71 


-0.83 


-0.90 


12.00 


0.29 


-0.23 


-0.51 


-0.73 


-0.81 


-0.92 


12.25 


0.07 


-0.28 


-0.51 


-0.79 


-1.50 


-2.00 


12.50 


-0.58 


-0.70 


-0.80 


-1.50 


-2.00 


-4.00 


12.75 


-0.78 


-0.70 


-1.50 


-2.00 


-4.00 


-9.00 


13.00 


-1.58 


-0.80 


-1.80 


-4.00 


-4.00 


-9.00 


13.25 


-4.00 


-4.00 


-4.00 


-9.00 


-9.00 


-9.00 



Table A2. Same as table IAT] but for model 4), and only showing 
cooling efficiencies for haloes more massive than W^^'^h~^MQ. 
For haloes less massive, the values arc the same as in table lAll 



Log(Mhalo) 


t = 0.80 


2.24 


3.38 


5.97 


10.27 


13.58 




2 = 7 


3.0 


2.0 


1.0 


0.3 





11.50 


0.71 


0.21 


0.02 


-0.32 


-0.41 


-0.67 


11.75 


0.43 


-0.15 


-0.45 


-0.71 


-1.83 


-1.90 


12.00 


0.29 


-0.23 


-0.51 


-0.83 


-2.51 


-2.52 


12.25 


0.07 


-0.28 


-0.70 


-1.00 


-4.00 


-4.00 


12.50 


-0.58 


-0.70 


-0.80 


-4.00 


-4.00 


-4.00 


12.75 


-9.00 


-9.00 


-9.00 


-9.00 


-9.00 


-9.00 


13.00 


-9.00 


-9.00 


-9.00 


-9.00 


-9.00 


-9.00 


13.25 


-9.00 


-9.00 


-9.00 


-9.00 


-9.00 


-9.00 



explicitly in Table \Kl\ The extrapolation and smoothing to 
the cooling efficiencies used in model of NW2010 can be 
seen by comparing Fig. lAll with Fig. 6 in NW2010, and 
also comparing the values listed in Table 1 with Table 6 
of NW2010. The cooling efficiencies shown in Fig. lAll and 
Table lAll are also applied to model 1), 2), 3) and 3b). 

Fig. IA2l shows the cooling efficiencies used for model 4) 
as presented in Sec. 2. 6, and table lA2l lists the explicit values 
of those efficiencies. 
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Figure A2. Cooling eificiency as a function of halo mass and 
cosmic time in model 4). The gray scale shows Log values of 
cooling efficiency in units of Log[Gyr~^]. X-axis is the cosmic 
time, with 13.7 Gyr corresponding to present day. 
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